Itinerant Ferromagnetism in a polarized two-component Fermi gas 
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We analyze when a repulsively-interacting two-component Fermi gas is thermodynamically stable 
in a phase separated state. We focus on the strongly polarised limit where the many-body problem 
can be described accurately in terms of repulsive polarons and the entropy of mixing. Phase diagrams 
as a function of polarisation, temperature, mass imbalance, and polaron energy as well as scattering 
length and range parameter are provided. When extrapolated to zero polarisation, our results are 
consistent with Monte-Carlo calculations. The lifetime of the polaron state is shown to increase 
significantly with the interaction range, which raises the prospects of realizing the elusive itinerant 
ferromagnetic phase with ultracold atoms at a moderately narrow Feshbach resonance. 



A mixture of two distinguishable species of fermions 
with strong repulsive interactions was predicted long ago 
by E. Stoner to undergo a transition towards an itiner- 
ant ferromagnetic state, where the two components phase 
separate to minimize the interaction energy [11. This is 
a strongly interacting many-body problem, notoriously 
difficult to solve, as such there are many competing and 
often contradictory theories. Moreover, its practical im- 
plementation is far from trivial and itinerant ferromag- 
netism has not yet been unambiguously observed. An 
exciting new setting to investigate this phenomenon is 
provided by ultracold atomic gases, since they offer great 
experimental flexibility and a precise determination of 
physical parameters [Jj. Though a recent experiment 
claimed the observation of itinerant ferromagnetism in 
a strongly repulsive Fermi gas [3], the current consensus 
is that fast decay precludes the realization of the ferro- 
magnetic state for atoms interacting via a broad Fesh- 
bach resonance [4]-[7]. However, a major step forward 
in the quest to produce long-lived repulsive Fermi gases 
was represented by the production of a ^Li-^'^K mixture 
at a moderately narrow Feshbach resonance, which leads 
to a much smaller decay rate [S]. This raises the hope 
that one can in fact observe itinerant ferromagnetism in 
atomic gases by choosing appropriate atoms and reso- 
nances. Previous theoretical studies of itinerant ferro- 
magnetism in the cold atom context employed mean-field, 
diagrammatic, and Monte-Carlo calculations [H [3 H)- 
[18] . The effects of a non-zero range of the interaction, 
which in light of their positive effect on the lifetime are 
highly relevant, were not examined in detail, and only in 
Refs. [3HS1II11 was the problem of decay considered. 

In this paper, we perform a thermodynamic analysis by 
comparing the free energies of two possible configurations 
of a two-component Fermi gas, the fully mixed phase, 
and the partially/fully separated (ferromagnetic) phase. 



A key point is that we focus on the limit of strong po- 
larisation where one can develop an accurate many-body 
theory in terms of repulsive polarons and the entropy of 
mixing. This allows us to calculate reliable phase dia- 
grams as a function of polarisation, temperature, mass- 
imbalance and polaron energy as well as scattering length 
and range parameter. The ferromagnetic region is shown 
to shift toward the BCS side with increasing range pa- 
rameter while the lifetime of polarons increases signifi- 
cantly. We finally perform a virial expansion, which con- 
firms that the critical temperature for ferromagnetism 
decreases with increasing range. The substantial lifetime 
improvement indicates that the itinerant ferromagnetic 
phase might finally be realised with cold atoms interact- 
ing via moderately narrow Feshbach resonances, which 
are readily available experimentally. 

We consider a two-component Fermi gas consisting of 
iVi atoms with mass mi and N2 atoms with mass TO2 
in a volume V . Two atoms of different species interact 
repulsively ("on the upper branch") via a short-range 
potential in the s-wave channel, whereas the intra-species 
interaction can be neglected. If the system is perfectly 
phase separated so that the 1-atoms occupy a volume Vi 
and the 2-atoms occupy a volume V2 — V — Vi with no 
mixing of the two species, the energy per particle is 

Escp = (1 - y)e,{N,/Vi,T) + ye^iN^/V^^T), (1) 

where y — N2/N is the fraction of the 2-atoms, and N — 
Ni + N2. The energy per particle of an ideal Fermi 
gas of Ncr (T-atoms in a volume reads 



where hix{z) is the polylogarithm. The fugacity ~ 
ex^{^ia /kBT), with the chemical potential, is de- 
termined by the density n^. = N„/Va via Ua = 
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-Li3/2(^^)/A^ where = {21^%^ /kBTrriaf/'^ is the 
thermal de Broghe wavelength. We recover the usual 
results = iEpal"^ + "^{keT / EFa?]/h for T < Tp^, 
and Ecr = 3kBT/2 for T ^ T^cr, where Epa- = kBTp„ = 
?i2(67r2n,)2/3/2m<,. 

Consider next a highly polarised homogeneous mixture 
with N2 ^ Ni. In this limit the system can be described 
as an ideal gas of 1-atoms, mixed with a gas of well- 
defined quasiparticles consisting of 2-atoms dressed by 
the 1-atoms. The quasiparticles are referred to as repul- 
sive polarons and they are the many-body equivalent of 
the repulsive branch of a Feshbach resonance. The energy 
per particle in this limit is accurately given by [13] 

£mix = {l~y)ei{N^/V,T)+ye2{N2/V,T)+v[l-vf'^E+. 

(3) 

Here Ej^ is the energy of a zero-momentum repulsive 
polaron in a Fermi sea of 1-atoms with density n = 
N/V and Fermi energy Ep = h'^kp/2mi, with kp = 
(GTT^n)^/^ [6]. In the low temperature regime, we may 
safely approximate E^- by its value at zero tempera- 
ture. The factor (1 — y)"^^^ is a rescaling of the po- 
laron energy taking into account the fact that the 2- 
atoms are immersed in a Fermi sea with lower density 
Ni/V = (1 - y)N/V. The effective mass of the re- 
pulsive polaron is close to TO2 [S], and we therefore take 
= ™2 for simplicity. Fixed-node Monte-Carlo cal- 
culations indicate that the expression Q is accurate for 
polarizations P = {Ni - N2)/{Ni + N2) > 0.5, corre- 
sponding to y < 0.25 US]- 

The thermodynamically unstable region of the homo- 
geneous mixed phase can be determined by applying the 
usual Maxwell construction to the free energy per particle 
/ = e—Ts, where s is the entropy per particle. Let us first 
focus on the case of equal masses mi = = m. The 
case of mi 7^ 7712 will be discussed later. For the phase 
separated state, equating the pressures of the two phases 
yields Ni/Vi = N2/V2 = n, the total density of the sys- 
tem. The energy per particle ([T| is then Sscp = e{n,T) 
independent of y. Here we have dropped the subscript 
a on Sa since mi = m2. We can therefore perform a 
Maxwell construction on the difference A/ = /,„ix — /sop 
which is given by 

A/ - {l-y)eini)+ye{n2)+yil-yf^''E+-e-TAs. (4) 

A major advantage of using the difference A/ is that 
we do not need to calculate Smix for the mixed state. 
Instead, we can use the entropy of mixing expression 
As = —ksiylny + {1 — y) ln(l — y)], which is the change 
in entropy corresponding to complete mixing of two sep- 
arated ideal gases. This expression reflects the purely 
combinatorial effects of mixing and we therefore expect 
it to be accurate for y ^ 0.25, where interactions between 
quasiparticles are negligible (ideal mixture assumption). 

In Fig. [TJ we plot the phase diagram obtained from a 
Maxwell construction using Q as a function of E+ and 
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FIG. 1. Phase diagrams in terms of the polaron energy E+, 
the polarization P, and the temperature T. The gas is mixed 
above the lines, and phase separated below. 

P for various temperatures (top), and as a function of 
P and T for various polaron energies E+ (bottom). We 
have only shown the diagram for P > 1/2 where our 
theory can be expected to be accurate. When T — 
and P — >■ 1, Fig. [T] shows that the system phase sep- 
arates when E+ > Ep. This simply reflects that the 
2-atoms cannot diffuse into a polaron state in the ideal 
gas of 1-atoms if the polaron energy is higher than the 
Fermi energy. For smaller polarisation P, phase sepa- 
ration occurs at a smaller polaron energy E^ since the 
system can separate into two partially polarised phases, 
which reduces the kinetic energy cost. Conversely, we 
see that phase separation is suppressed at higher tem- 
peratures due to the entropy of mixing. Note that the 
phase diagram in Fig. [l] is generic in the sense that it 
is based only on the existence of well-defined repulsive 
quasi-particles with energy which has been verified 
experimentally [S] , and on the ideal mixture assumption. 

To make a direct connection with experiments, we now 
calculate E+ in terms of the scattering parameters of the 
interaction: the s-wave scattering length a and the range 
parameter R* . At the many-body level, a small/large 
kpR* , corresponds to a wide/narrow Feshbach resonance. 
The range parameter is R* — / {2mrahgS fiAB) > 0, 
where Cbg is the background scattering length, is 
the differential magnetic moment, and AB the mag- 
netic width of the Feshbach resonance. A very attrac- 
tive feature of the high polarization limit is that it is 
possible to calculate E+ precisely even for strong inter- 
actions [51 |S1 [TM^ . The details of the many-body cal- 
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FIG. 2. Phase diagram for equal masses, at a broad reso- 
nance with kpR* = (thin lines), and at a narrow one with 
kpR* ~ 1 (thick). The gas is mixed above the lines, and 
phase separated below. The black markers in the graph indi- 
cate the IFM transition at P = T = i?* = as found using 
(from left to right) Monte-Carlo calculations [11H13) . second 
order perturbation theory [9], and mean field (Stoner) the- 
ory d]: = {0.9, 1.1, 7r/2}. 



culation of as a function of a and R* , which includes 
one particle-hole excitations from the Fermi sea of the 
majority atoms, can be found elsewhere [6, 8^, ^ . 

The resulting phase diagram in terms of the scattering 
parameters is shown in Fig. [2] We consider two different 
values of the range: kpR* = (thin lines) describing a 
wide Feshbach resonance, and kpR* — 1 (thick lines) cor- 
responding to the intermediate value of the experiment in 
Ref. ^ . We observe that a non-zero effective range shifts 
the phase separated region toward the BCS-regime (the 
right side of the plot), consistent with the shift of the po- 
laron/molecule crossing reported in Ref. 8J. As in Fig. [T] 
we have drawn the phase boundary lines only within the 
regime of validity of the polaron theory, i.e., P > 1/2. 
With increasing —{kpa)~^, we have furthermore termi- 
nated the lines where the polaron Ansatz fails as the 
polaron decay rate Tpp becomes too large, i.e., when 
Tpp/Ep > 0.25. The issue of the polaron decay will 
be discussed in more detail below. For comparison, in 
Fig. [2] we also give the critical values for the ferromag- 
netic transition for the balanced case P = calculated 
by various zero-temperature theories pQ IHl fTTHTB] . Due 
to the symmetry of the phase diagram for equal masses, 
the phase boundary must cross the P — axis vertically. 
As can be seen from Fig. [2j this significantly restricts the 
range of possible extrapolations of our theory from its 
range of validity to smaller values of \P\. In particular, 
the extrapolation of our theory to P = is consistent 
with the r = critical coupling strengths predicted by 
the Monte-Carlo calculations [ll,-,-13j which indicates the 
accuracy of our approach. Conversely, it is difficult to 
reconcile with mean-field Stoner theory [T]. 

Since a long lived repulsive polaron state has been re- 



cently achieved for a mixture of ^Li and ''"K atoms [S], we 
turn now to study the effects of unequal masses toi 7^ 7712 . 
For simplicity, we consider T = and the limit P ^ 1 
where it is sufficient to compare the energy of the mixed 
phase with that of the perfectly phase separated state, 
given respectively by (|3| and ([T]). For the perfectly 
phase separated state, the pressures Pa of the two phases 
must be equal; using a [Nfj /Va)^^^ /m„ for an ideal 
Fermi gas, we find Ni/Vi = n[y(mi/m2)^/^ + 1 — y] and 
N2/V2 = n[{\ ~ y){m2/miY/^ + y\. To first order of y, 
Eq. ([I]) becomes 



Escp = Epi{n) 3/5 + y(mi/TO2)^/^ 



(5) 



which is accurate when y <C 3|(mi/m2)'^/^ — To 
the same order of y, Eq. ([3| is of the form Emix = 
Epi{n) [3/5 -I- y{E+/Epi{n) — 1)]. The homogeneous 
mixture is unstable when Smix > £,s 



E+ > 



nil 

1112 



3/5 



Epi{n) = 



scp 

which gives [6] 
(6^n)2/3 



o 3/5 2/5 ■ 

2m2 m-^ 



(6) 



This result demonstrates that phase separation is favored 
by making one of the atom species very heavy, since the 
kinetic energy cost then decreases. This effect was also 
discussed in Refs. [15j [18] . 

Previous studies showed that strongly repulsive Fermi 
gases prepared close to a broad Feshbach resonance suf- 
fer from severe decay toward the "lower branch" [5HZ] ■ It 
is therefore important to address the question: are there 
conditions under which a strongly repulsive Fermi gas has 
a sufficiently long lifetime for the observation of the fer- 
romagnetic state? To examine this, we plot in Fig. [3j^a) 
the critical value of the interaction parameter 1/kpa for 
phase separation at T = and P — 1 for different mass 
ratios obtained from (|6]), and in Figure |3]jb) we plot the 
corresponding two-body decay rate Fpp of the repulsive 
polaron at the critical coupling strength. The rate Tpp 
is calculated using a diagrammatic method including ef- 
fects from a non-zero range R* , whose results are in good 
agreement with experimental data [5J |H] • Figure [s] shows 
an important result: a narrow resonance with large kpR* 
gives rise to a longer polaron lifetime than a broad one. 
For instance, for a mixture of a few '^^K atoms in a gas 
of ^Li atoms, the mass ratio is m2/mi — 40/6, and the 
polaron lifetime increases by a factor ^ 10 at the criti- 
cal coupling strength for phase separation if the range of 
the atom-atom interaction is kpR* = 1 instead of zero. 
Furthermore, we find that a large mass ratio 7712/ mi de- 
creases the decay rate significantly compared to the case 
of equal masses. These results demonstrate that the de- 
cay of the repulsive phase towards the attractive branch 
may be strongly reduced, and as such the prospects of ob- 
serving itinerant ferromagnetism increased, by employing 
heavy impurities in a bath of light atoms at a moderately 
narrow Feshbach resonance. Note that our results indi- 
cate that phase separation of an initially mixed phase is 
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FIG. 3. (a) The critical coupling strength for phase separation 
for T — and P = 1 as a function of R* . (b) The decay 
rate Fpp of the repulsive polaron at the transition. Blue: 
milm\ = 1, black: 40/6, and red: 6/40. 




unlikely due to fast decay (larger than typical trap fre- 
quencies), but that an initially separated phase should 
be long-lived. To investigate further the stability of the 
phase separated state, one could calculate the tunneling 
rate of atoms between the two phases taking into account 
not only the lifetime of the repulsive polaron but also its 
quasiparticle residue and the molecule-hole continuum. 

We finally examine the high temperature regime, 
where one can obtain reliable results using the virial ex- 
pansion. To second order, this expansion gives for the 
partition function of the two-component Fermi gas 



Z = 1 




(7) 



with M — nil + 7712, l/^r = + l/"^2j — 

{2TTn^/2kBTmr)^/^, and Am = {2TTh^ /keTM)^/^ . When 
only the s-wave "upper branch" of excitations is taken 
into account, the second virial coefficient 62 is [23] 



bo = 



1 



nirT 



dk 



kSs{k)t 



(8) 



The s-wave phase shift is given by coti5s(fc) = 
— 1/ka — R*k. At a broad resonance one has &2 — 
(X^rli^af \^ ^ erf (Ar/2V7Fa) - 261(a)] /2, whose mini- 
mum value is —1/2 at 1/a = 0+. For a very narrow 
resonance, we find 62 ~ —Xr/^TiR* at the unitary point. 

The free energy per particle of the mixture /juix = 
-ksTlnZ/N + (1 - v)iii + is 

/mix = -fcsT {1 - (1 - y) ln[(l - y)nXl] - y HynXl) 
+y{l - y)nXlb2} , (9) 

where we have used Zo- = A^(no- — A^nin2&2) obtained 
from ([7]) . The phase diagram will have a phase separated 
region when 

a^/mix/fesT = (1 - y)-^ + y-' + 2nX% < 0, (10) 

for some values of y, which is possible only if 62 < 0. For 
equal masses, the region of phase separation derived from 
the Maxwell construction is bounded by the concentra- 
tions which minimise /mixj i-e., 



FIG. 4. Phase diagram at the unitarity point = for 
equal masses as obtained from the virial expansion to second 
order: kpR* = (thin line) and kpR* = 1 (thick). The gas 
is mixed above the lines, and phase separated below. 



The phase diagram obtained from (11) is plotted in 
Fig. [4] for (kpa)^^ = 0, toi = m2, and kpR* — 
or kpR* — 1. Since there is no kinetic energy cost of 
phase separation for high temperatures, this phase di- 
agram is a result of the competition between the re- 
pulsive interaction energy and the entropy of mixing. 
The predicted critical temperature for phase separation 
unfortunately is too low for the virial expansion to be 
reliable. Indeed, the largest transition temperature is 
jimax 0.66Tp obtained at unitarity for a broad reso- 
nance where 62 — —1/2. The calculation nevertheless 
indicates that a large range decreases the critical tem- 
perature for phase separation, consistently with the low 
temperature polaron calculation. 

To conclude, we derived detailed phase diagrams of a 
two component Fermi gas in the limit of strong polari- 
sation, where the effects of the interactions can be accu- 
rately described in terms of the repulsive polaron energy 
and the entropy of mixing. This phase diagram was then 
expressed in terms of the scattering length and range pa- 
rameter using a many-body theory known to be reliable 
in the strongly polarised limit. The ferromagnetic region 
was shown to move towards the BCS side with increas- 
ing range of the interaction. Importantly, a large range 
furthermore increases substantially the polaron lifetime, 
raising hopes for the realisation of itinerant ferromag- 
netism in a physical system. 
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